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In a recent paper [Ferrari et ai, Phys. Rev. E 69, 035102(R) (2004)], the scaling law of the 
fluctuations of the step limiting a crystal facet has been computed as a function of the facet size. 
Ferrari et al. use rigorous, but physically rather obscure, arguments. Approaching the problem from 
a different perspective, we rederive more transparently the scaling behavior of facet edge fluctuations 
as a function of time. Such behavior can be scrutinized with STM experiments and with numerical 
simulations. 
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In recent years it has become possible to probe quanti- 
tatively with STM the detailed atomic-scale fluctuations 
of steps near facet edges, most notably illustrated in ex- 
tensive, painstaking measurements near (111) facets of 
Pb microcrystallites on a Ru support 0, Bonzel's 
review provides a thorough and lucid account of these 
investigations [3j. Accordingly, it is an opportune time 
for theoretical examination of such systems. 

In a very recent intriguing paper Q , Ferrari et al. have 
computed the scaling of equilibrium fluctuations of an 
atomic ledge bordering a crystalline facet. These authors 
find that the step-edge width w scales as w ^ L 1 / 3 with 
the linear size L of the facet. This result differs from what 
is expected, and actually found, for the step bordering a 
2D island, which performs a random walk so that w ~ 
L 1 / 2 . Ferrari et al. claim that the origin of the unusual 
L 1 / 3 scaling lies in the step-step interactions between the 
facet ledge and the neighboring steps. 

Ferrari et aVs calculation is based on the use of free 
fcrmions, transfer matrix, random matrix properties, and 
specific models; it does not address the question of the 
time behavior of step fluctuations. 

In the present Communication, we approach the prob- 
lem from the perspective of a continuum-equation de- 
scription of a faceted crystal. The most easily accessible 
experimental quantity is the step autocorrelation func- 
tion G(t) = ([x(t) — x(Q)] 2 ), which is expected to have a 
power-law behavior at short times: G(t) ~ t 2 ^. Hence, 
we compute the scaling of the ledge fluctuations with 
time, using both the continuum-equation approach and 
the simple arguments developed in Pimpinelli et al. [jj. 

Since the development in Ferrari et al. Q is rather ob- 
scure physically obscure, we readdress the problem from 
a different point of view, based on continuum equations 
and scaling. Our approach shows that the mean square 
width of a fluctuating crystal surface next to a facet scales 
as ((8z) 2 ) ~ £ 2 / 3 , I being the length scale in the radial di- 
rection. We obtain then the time behavior of the surface 



fluctuations, that will then be compared to that obtained 
from a different, more qualitative approach (see below). 

The projected free energy of a surface near a facet be- 
low its roughening temperature is given by 0, 0, 

f(<f>) = 7b CO + ^jp- tan + g(9, T) tan 3 <f>, (1) 

where h is the step height, 7o(? 1 ) is the terrace contribu- 
tion, (3(0, T) is the step free energy per length, and g the 
"step interaction parameter" • The angles <j> and are 
the angle of the surface relative to the facet and of the 
steps relative to an arbitrary direction, respectively. 

The chemical potential for a crystallite in cylindrical 
coordinates (r, 0, z) reads 

= v0(T) | 3vg(T) fdz{r)\ 2 d (dz{r)\ 2 

hr r \ dr J dr \ dr J 

(2) 

where the 6 dependence of the step free energy and of 
the interactions has been neglected for simplicity. Here v 
is the atomic volume, with v/h then the atomic surface 
area. It is straightforward to verify that the Pokrovsky- 
Talapov equilibrium shape 0] is recovered from Eq. J2J. 
In the present geometry, the equilibrium shape is a stack- 
ing of circular layers separated by steps. 

Consider the top step of the stack. Locally the ledge 
performs a random walk, so that each time that the ledge 
moves forward or backward by one unit, r — po±Sr, the 
local surface height increases or decreases by one unit, 
respectively, z{r) = z$ ± 8z. (Note that the variable Sz 
is the continuum translation of h.) In Eq. J2J the first 
term is the Gibbs-Thomson contribution coming from the 
curvature of the layer, the second term represents the 
variation in the interaction energy due to the change of 
the length of the ledge when a single atom is removed or 
added. The last term represents the change of interaction 
energy when the step-step distance in the curved part of 
the crystal is varied. 
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Thus, the contribution which dominates in Eq. for 
small fluctuations of the surface is the second one, since 
the curvature and the step-step distance are weakly af- 
fected by local step motions. Both the first and the last 
term in Eq. (J2J can then be neglected. After next setting 
1/r in the second term in Eq. @ equal to l/p® to first 
order and letting d(5z)/dt = —K\I 2 \i + i] c (where K is 
a transport coefficient, and -q c is a stochastic term de- 
scribing mass-conserving random atomic motions at the 
interface), we arrive at the equation 



dt 



•v.'/< :Q /£ 



Pa 



dt 



+ r)c(r,t). (3) 



Eq. © looks like the nonlinear part of the so-called 
"Montreal" or conserved KPZ model [ljj. However, in 
the original model the nonlinearity is absent in equilib- 
rium, while here it stems from the contribution of step- 
step interactions to the equilibrium chemical potential. 

Since Eq. is nonlinear, it cannot be solved ana- 
lytically. Hence, we limit ourselves to a power-counting 
evaluation of the scaling exponents, following Hentschcl 
and Family ^lj. Taking a characteristic length scale £ 
as the unit, we can estimate the root- mean-square value 
of the time derivative of the fluctuations as Szg/ti, and 
of their spatial derivative as Sz 2 /£ . The root-mean- 
square value of the conserved noise term is estimated as 
r] c ~ l/(e 2 S e t e ) 1 / 2 , where Si = y/P + z'j is the length 
of the fluctuating interface. For small-amplitude fluctua- 
tions, as is the case here, Se ~ £, so that r\ c ~ l/(£ 3 ^) 1 / 2 

Equating the time derivative to the noise term yields 



Sze 



(4) 



Equating the spatial derivative to the noise term yields 
Sz^^/t 1 / 2 . (5) 
From Eqs. 10} and (J5J we finally get 

Szf ~ £, (6) 

which has the same form as Ferrari et a/.'s result. Note 
that the dynamics of step fluctuations is also affected: 
From Eqs. (0J and JJJ, e.g., we obtain 



x A/ 11 
bzi ~ t/ , 



as well as 



£~t 



3/11 



(7) 



(8) 



Our argument allows extensions to nonconservative 
ledge fluctuations (i.e. fluctuations driven by attachment- 
detachment of atoms to and from the ledge). Indeed, 



if ledge fluctuations are driven by atom attachement- 
detachment, with kinetic coefficient K, Eq. J2J) has to 
be replaced by 



dt 



Pa 



dt 



Proceeding as before, we can estimate the time deriva- 
tive of the fluctuations as Szi/tg, their spatial deriva- 
tive as SzjJJ 2 , and the (nonconserved) noise term as 
1 / '{ItiY I 2 11]. In the latter relation we assume again 
small-amplitude ledge fluctuations, consistently with the 
hindering action of neighboring steps. 

Equating the time derivative to the noise term yields 

sz e ~ 4 /2 /^ 1/2 - ( 10 ) 

Equating the spatial derivative to the noise term yields 
Szt-iW/t 1 / 2 . (11) 
From Eqs. lfTT)jl and fTTJ) we obtain 

u ~ £ 5/3 , (12) 

as well as 

Sz e ~ e 1 / 3 . (13) 
Together Eqs. (JHJ and O yield 

Sze ~ t 1 / 5 , (14) 

as found above. 

Note the Eq. looks like the KPZ equation with- 
out the linear term. However, the nonlinearity stems 
here from equilibrium, small-amplitude fluctuations of 
the step edge, and the resulting dynamics is different 
from that of the KPZ model. The latter may be recovered 
from Eq. @ by assuming large-amplitude fluctuations, so 
that Si ~ Szi, and rji ~ l/(Szitg) 1 ^ 2 . 

Equating the time derivative to the noise term yields 



X 3 / 2 +V2 

Sz/ ~ t/ . 



(15) 



Equating the spatial derivative to the noise term yields 
Sz 5 / 2 ~ £ 2 /t] /2 . (16) 
From Eqs. I|15|l and l|16f) we obtain 

ti ~ £ 3/2 , (17) 

as well as 

Szi ~ £^ 2 , (18) 

reproducing the power laws characteristic of the KPZ 
model. Note that in the limit of unhindered, large am- 
plitude fluctuations, the random walk scaling w ~ £ 1/2 is 
recovered, as expected. 
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Let us assume now that the results of Ferrari et al. 
hold, and let us use Pimpinelli et aVs argument |5j for 
finding the time scaling. The argument computes the 
width 5x of a fluctuation of linear size £ is found by letting 
the area of the fluctuation, £Sx, equal to the fluctuation 
8N of the number N(t) of particles involved in the mass- 
transport process responsible for the fluctuation, during 
a time t. In turn, N(t) is estimated as 



N(t) 



( <\ 



lL s t, 



(19) 



where c eq /r* is the rate of the transport process, c eq be- 
ing the equilibrium particle density, and £L S the surface 
area feeding the fluctuation ||. The assumption that 
mass transport is conservative, which is consistent with 
Eq. (3) above, yields c eq /r* c eq D e /£ 2 , where D e is the 
edge diffusion coefficient. Also, £L S » al, a being of the 
order of the lattice spacing. Thus, 

N(t) « tc eq D e a/£. (20) 

Letting SN = y/N, we now find 

{5x) 2 £ 2 w (SN) 2 w tc eq D e a/£. (21) 

Assuming the scaling relation w ~ I 1 / 3 yields finally 

* ~ £ 11/3 , (22) 



6x(t) « t 



1/11 



(23) 



so that G(t) ~ t 2 / 11 . 

As in Ref. 0, we can also compute the time behavior 
of a fluctuation with non-conserved kinetics. This can 
be done by letting c eq /r* w kc eq , with k an appropriate 
kinetic coefficient for atom detachment / attachment from 
and to the step edge. Then, using Eq. (|19l) with L s sa a 
yields 



(Sx) 2 £ 2 w fcc e ^crf. 



(24) 



Again letting 5a; ~ i 11 / 3 yields £ 5 ^ 3 ~ 4, eventually 

&c - t 1 / 5 . (25) 



Thus, the temporal scaling laws turn out to be identical 
to those computed above for the fluctuations of the facet 
edge in the continuum equation approach. 



In summary we have shown how the powerful general 
arguments used a decade ago to launch the systematic ex- 
ploration of fluctuations of steps on vicinal surfaces can 
be extended to examine the fluctuations of steps near a 
crystalline facet. We find rich and varied behavior that 
we hope will stimulate closer examination of fluctuation 
phenomena near step edges by both experiment and sim- 
ulation. 
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